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Abstract 

The pp ppr]' reaction is investigated within a relativistic meson-exchange 
model of hadronic interactions. We explore the role of nucleonic and mesonic, 
as well as the N* resonance currents, in producing r/' mesons. In order to learn 
more about the production mechanisms, new measurements in the energy 
region far from the threshold are required. 
PACS: 13.60Le, 25.10+s, 25.40-h 
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I. INTRODUCTION 



With the advent of particle accelerators in the few GeV energy region, heavy meson 
production in hadronic collisions has attracted increasing attention in the past few years. 
In particular, heavy meson production in nucleon-nucleon (NN) collisions at near-threshold 
energies is of special interest, not only because it is suited for extracting information on a 
few lowest-order multipole amplitudes, but also because it is considered to provide impor- 
tant information on the short distance behavior of the A^A^ interaction. Due to the large 
momentum transfer between the initial and final nucleons, these reactions at near-threshold 
energies necessarily probe the NN interaction at short distances. In the present work we 
concentrate on the t]' meson production in proton-proton (pp) collisions. 

Among heavy mesons, the rj' meson is of particular interest for various reasons. The r]' 
meson is thought to couple strongly to gluons via the QCD anomaly coupling r]' ^ g + g 
imi. Also, an unexpectedly large branching ratio measured recently for the inclusive decay 
of beauty particles, B ^ t]' + X has been interpreted as possible evidence for the strong 
coupling of rj' meson to the gluonic components [Q. It would then be conceivable that the 
pp —>■ pprj' reaction might probe the gluon content of the f]' meson via its coupling to gluons 
emitted from the quarks exchanged between two interacting nucleons. This mechanism 
would be complementary to the vector meson-exchange current mechanism for producing rj' 
mesons. 

One of the properties of the t]' meson of extreme importance is its yet-poorly-known 
coupling strength to the nucleon. This has attracted much attention in connection with 
the so-called " nucleon-spin crisis" in polarized deep inelastic lepton scattering 0]. The 
NNt]' coupling constant, gNN-q', may be related (through the axial vector coupling using the 
Goldberger-Treiman relation) to the quark helicity contribution to the spin of the proton . 
Therefore one can argue that gNNri' tells us about the total spin of the nucleon carried by 
its constituents; conversely, the quark contribution to the spin of the nucleon would tell us 
about gNNr]'- So far there is no direct experimental measurement of gNNri'- Recently, cross 
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sections for the pp pprj' reaction near threshold have been measured by the COSY-11 
and SPESIII collaborations. This reaction may offer an opportunity to determine 
this coupling from a direct emission process of r]' by a proton. Of course, other production 
mechanisms, such as meson exchange and nucleon resonance currents, must be taken into 
account before a quantitative determination of gNNri' is possible. The major aim of the 
present work is to explore the roles of various production mechanisms. 

The theory of rj' production in pp collisions is still in its early stage of development 
UsHTTj. In this work we investigate this reaction using a relativistic meson-exchange model 



of hadronic interactions. In section 2 we outline the formalism for calculating the production 
amplitude. The final state interaction is known to play a crucial role in the production of 
particles near threshold energies in A^A^ collisions ; therefore our formalism includes the 



pp final state interaction explicitly. The Coulomb correction in the pp final state interaction 
is also known to be important when these protons have small energies [p!3| , p!4| . In the present 



work it is treated exactly and reduces the calculated cross section by as much as a factor 
of two for energies close to threshold. The initial state interaction is taken into account 
via a reduction factor determined from the available phase shifts and inelasticities in 



the threshold incident energy region. In section 3 our rj' meson production currents are 
constructed. We consider the nucleonic, mesonic and resonance currents. The roles of these 
currents are explored in section 4. Our results are summarized in Section 5. 

II. FORMALISM 

The formalism used in the present work is essentially the same as that employed in Refs. 



T7| , p!8[| for studying the production of vector mesons. It is based on a relativistic meson- 
exchange model of hadronic interactions in which the transition amplitude is calculated in 
the Distorted Wave Born Approximation, taking explicitly into account effects of the final 
state interaction. We write the transition amplitude describing the p+p —>■ p+p + rj' process 
as 
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M=<0/|(T)-)^zG/ + l)J|0i> , (1) 

where 0jj denotes the four-component unperturbed pp wave function in the initial (i) and 
final (/) state. Tj"-* is the final state pp T-matrix. Gf stands for the two-nucleon propagator 
and J is the ?7'-emission current defined in the next section. 

The T-matrix used in our calculation is generated by solving a three-dimensional reduced 
Bethe-Salpeter equation (the Blankenbecler-Sugar equation) for a relativistic one-boson- 
exchange A^A^ potential V, i.e., 

T = V + VtGBBsT , (2) 

where Gbbs denotes the Blankenbecler-Sugar (BBS) two-nucleon propagator. In this work 
we employ a slightly modified version [|T^] of the Bonn B NN model as defined in Table 
A.l of Ref. for constructing the potential V. This modification has been made in order 
to reproduce the pp low-energy parameters rather than the pn low-energy parameters while 



preserving other features of the Bonn B interaction model for describing pp scattering |]I9 



It should also be mentioned that each nucleon-nucleon-meson (NNM) vertex in the A^A^ 



potential is modified by a form factor of either monopole or dipole form. We refer to ||20 
for further details. Furthermore we note that the two-nucleon propagator Gf appearing in 
Eq.(P is, for consistency, also chosen to be the BBS propagator, Gf = Gf-BBS- 

In Eq.(|l]) the initial state interaction is neglected. For heavy meson production such as 
1]', its effect on the production cross sections near threshold may be taken into account by 
multiplying the calculated cross sections using Eq.(|l]) by a factor |T5 | 

K = VaCOS^iSa) + - r]a]^ (3) 

where a stands for the quantum numbers specifying the corresponding initial A^A^ state and 
6a and rja denote the corresponding A^A^ phase-shift and inelasticity, respectively, at the 
nucleon incident energy. Although this reduction factor is only meant to account for the 
gross effect of the initial state interaction, it works well when confronted with a calculation 
where the initial state interaction has been treated explicitly In the near threshold 



energy region the relevant initial pp state is a = ^Pq. For rj' production the threshold 
incident energy is about Tiab = 2.404 GeV. Using the values 5(3pg) = —50° and ?7(3pg) = 0.75 
from a partial wave analysis ||22|, we obtain the reduction factor of = 0.33 for the cross 
section. An alternative way to account for effects of the initial state interaction is to absorb 
these effects in the form factors at the t]' meson production vertices as has been done in 
Refs. OJl^. We choose the former method in the present work. 



III. PRODUCTION CURRENTS 

Within our hadronic model of strong interactions, the T^'-emission current J in Eq.(|I]) 
consists of a sum of the baryonic and mesonic currents. The baryonic current is further 
divided into the nucleonic and nucleon resonance (A^*) currents, so that the total current is 
given by 

J Jnuc ~l~ Jres ~l~ Jmec • (4) 

The individual currents are illustrated diagrammatically in Fig. |I|. In the following subsec- 
tions we construct each of these currents. 



A. The Nucleonic Current 

The nucleonic current is defined as 

Jnuc= Yl (Xj^SjU + UiSjTj) , (5) 
i=i,2 

with Tj denoting the NNrj' vertex and 5*^ the nucleon (Feynman) propagator for nucleon j. 
The summation runs over the two interacting nucleons, 1 and 2. U stands for the meson- 
exchange A^A^ potential. It is, in principle, identical to the potential V appearing in the 
A^A^ scattering equation, except that here meson retardation effects (which are neglected in 
the potential entering in Eq.(|^)) are kept as given by the Feynman prescription. 

The structure of the NNrj' vertex, F^-, in Eq.(|^) is derived from the Lagrangian density 
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a + ^-:^l^d\ r,\x) 



2m 



N 



(6) 



where gNNri' denotes the NNt]' couphng constant and A is the parameter controlhng the 
pseudoscalar(ps) - pseudovector(pv) admixture. ?7'(x) and \1/Ar(x) stand for the t]' and nu- 
cleon field, respectively; ttin denotes the nucleon mass. 

As mentioned in the introduction, the coupling constants QNNri' and the ps-pv mixing 
parameter A are poorly known at present. The predictions for gNNrj' range anywhere from 
|]; an estimate based on the dispersion method even gives gNNrj' consistent 



1.9 to 7.5 



with zero p5|. Zhang et al. ||2^ in their analysis of the photoproduction on protons of rj' 
mesons used A = 1 for the ps-pv mixing parameter. Bernard et al. 0, in their analysis of 
the pp —>■ pprj' process, extracted a value of A = 0.4 ± 0.1 in conjunction with the value of 
dNNrj' = 2.5 ± 0.7 determined from a recent measurements of AS in deep inelastic lepton- 
nucleon scattering p6[. The latter quantity is related to qnnt^' 0- SU(3) symmetry, together 
with the OZI rule [^, relates the NNrj' coupling to the NNr] coupling: gNNrj' = dNNri 
tan(ap), where ap = 9p — 6p(^ideai) — —45° denotes the deviation from the pseudoscalar 



ideal mixing angle. With the value of gNNr] = 6.14 used in A^A^ scattering analysis [pO[ 



we 



then have gNNr^' — 6.1 which is close to the upper end of the predicted range mentioned 
above. The value of gNN-q = 6.14 together with the rj — rj' mixing angle of 6p ~ —9.7°, as 
suggested by the quadratic mass formula, and the NNtt coupling constant of g^Nj^ = 13.45 
leads to the ratio D/F ~ 1.43. This is not too far from the value of D/F = 1.73 extracted 
from a systematic analysis of semileptonic hyperon decays p8|. In the present work we use 
the value of gNNrj' = 6.1. We shall consider both of the extreme values of the parameter A 
in Eq.(^, i.e., A = and 1. Note that rj' is not a Goldstone boson. Consequently there is, 
a priori, no constraint on the ps-pv admixture. 

The NNt]' vertex derived from Eq.(^) should be provided with an off-shell form factor. 
Following Ref. |T^, it is assumed to be of the form 



where denotes the four-momentum squared of either the incoming or outgoing off-shell 
nucleon. We also introduce the form factor given by Eq.(0) at those NNM vertices appearing 
next to the r/'-production vertex, where the (intermediate) nucleon and the exchanged mesons 
are off their mass shell (see Fig. |I|). Therefore, the corresponding form factors are given by 
the product F^{p)FM{q%i), where M stands for each of the exchanged mesons between 
the two interacting nucleons. The form factor FM{q\.i) accounts for the off-shellness of the 
exchanged meson and is taken consistently with the NN potential used for generating the 
T-matrix. In Ref. |ll8|] cutoff parameters in the range of A^v = 1-17 — 1.55 GeV have been 
used in the description of 0- and cj-meson productions in pp collisions. In order to reduce 
the number of free parameters, here, we adopt the value of A^r = 1.2 GeV , which is also the 
value adopted at the NN*rj' vertex in the construction of the resonance current in the next 
subsection. 



B. The Resonance Current 

Currently there are no well established nucleon resonances that decay into Nrj' . In 
Ref. the strong peak in the total cross section of the rj' photoproduction off protons 
close to threshold has been attributed to the Di3(2080) resonance. On the other hand, 
a recent multipole analysis of the r]' photoproduction [1^ indicates that this reaction is 
dominated by an 5*11 and a Pu resonance. In Ref. |T^, apart from the resonance dominance 
assumption — which probably leads to an overestimation of the total decay width of the 
resonances considered, it is also assumed that these resonances decay only into prj' and pvr". 
The latter decay channel also effectively accounts for all other decay channels that were 
not considered explicitly in the analysis due to lack of information. Although this may 
not affect appreciably the description of the rj' photoproduction reaction, the pp — > ppr]' 
process can be very sensitive to such an assumption. Since the pion couples strongly to the 
nucleon, an overestimation of the pvr" partial decay width is likely to lead to a considerable 
overestimation of the cross section due to these resonances in the pp — > ppr]' reaction. In 
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addition to all these uncertainties in the extraction of the resonance parameters, there is also 
the possibility of a threshold cusp effect, as discussed in Ref. that might explain the 
observed behavior of the photoproduction cross section close to threshold in the absence of 
any resonance. All these issues require further careful considerations. With this background 
in mind we will explore the role of the 5'ii(1897) and Pii(1986) resonances as determined in 



Ref. \16 



The resonance current, in analogy to the nucleonic current, is written as 



Jres — X! [rjN*r]''iSN*UN* + UN*iSN*'rjN*r]'J ■ (8) 

j=l,2 N* 

Here TjN*r)' stands for the NN*ri' vertex function involving the nucleon j. Sn*{p) = + 
m7v)/(P^ — m^, + imiy*T^*) is the N* resonance propagator, with rriiy* and T^* denoting 
the mass and width of the resonance, respectively. The summation runs over the two 
interacting nucleons, j = 1 and 2, and also over the resonances considered, i.e., A^* = 
5ii(1897), Pii(1986). In the above equation Un* (Un-) stands for the A^A^ A^A^* {NN* 
NN) meson-exchange transition potential. It is given by 

Un* = rAriv*^(g^)iA^(g^)rAfAr^(g^) , (9) 

where A^(g^) denotes the (Feynman) propagator of the exchanged pion with four-momentum 
q-jt and Tj^iNT^lq^^) the NNti vertex. The latter is taken consistently with the A^A^ potential V 
appearing in Eq.(|D. We note that there is also an additional contribution to Un* from the t]' 
exchange in Eq.(P). However this is negligible compared to the pion-exchange contribution. 

Following Ref. |3^, the NN*r]' and NN*7t vertices in Eqs.(P,|^) are obtained from the 
interaction Lagrangian densities 

^Nl*r,'i^) = -9nnw^n{x) I [^AFW + f ^"^ FW^a^I r]'{x)] ^n*{x) + h.c. 
' [[ \mN*±mN) ^ J J 

1^^il*A^) = -gNN*.^N{x) I Lf(±) + ( ^"^ \ F^S'^l r ■ 7f(x)| ^N'{x) + h.c. , (10) 



where 7f(a;) and \1/Ar. (x) denote the pion and nucleon resonance field, respectively. The 
upper and lower signs refer to the even(+) and odd(— ) parity resonance, respectively. The 
operators F*^^^ and F|^^^ in the above equation are given by 



7/^75 



r(-^ = i,r(-) = 7.. (11) 

The parameter A in Eq.(^) controls the admixture of the two types of couphngs: ps (A = 1) 
and pv (A = 0) in the case of an even parity resonance and, scalar (A = 1) and vector 
(A = 0) in the case of an odd parity resonance. On-shell, both choices of the parameter A 
are equivalent. Hereafter the mixing parameter A for the NN*'n' vertex is fixed to be A = 0. 
We consider the choices A = and 1 in the NN*!]' vertex, however. 

The coupling constants Qnn't]' and Qmn'-k are determined from the extracted decay widths 



(and masses) of the resonances from the r]' photoproduction on protons |T6 



{ms,, , Ts,,) = (1.897 ±0.050^°:°;^° , 0.396 ± 0.155+°;°^^) GeV 

{mp,, , Tp,,) = (1.986 ± 0.026t°:°^° , 0.296 ± 0.100+°;°?°) GeV , (12) 



with partial decay widths 



r(5'n ^ Pv') = 0.0515,, , r(5n ^ pvr°) = 0.95rs,, 

r(Pn ^ Pv') = 0.2515,, , r(Pn ^ pvr") = 0.75r5„ . (13) 
Using these values we obtain 

QNSiiri' = 2.9 , QNSiin = 2.4 
gNPiiv' = 11-7 , gNPim = 5.2 . (14) 

Here we assume the coupling constants to be positive. We mention in advance that the 
interference of the resonance current with the rest of the currents (nucleonic and mesonic 
currents) is practically negligible. 

In complete analogy to the nucleonic current, we introduce the off-shell form factors at 
each vertex involved in the resonance current. We adopt the same form factor given by 
Eq.(|^, with mjv replaced by m^* at the NN*t]' and NN*it vertices, in order to account for 



the off-shellness of the N* resonances. We note that Zhang et al. have also employed 



this form for the form factor at the NN*r]' vertex in their study of t]' photoproduction. In 
order to account for the off-shellness of the exchanged pion (see Eq.(y)), the NN*'n' vertex 
is also muhiphed by an extra form factor of F-^^ql) = (A^ — m^)/(A^ — ql), the same form 
used at the NNn vertex in the construction of the A^A^ potential V entering in Eq.(^), with 
the cutoff parameter of A^ = 0.8 GeV. We also use the same form factor at the NNn vertex 
in Eq.(l). 

C. The Mesonic Current 

For the meson-exchange current we consider the contribution from the vvrj' vertex with 
V denoting either a. p or uj meson. As we shall show later, this gives rise to the dominant 
meson-exchange current. The vvrj' vertex required for constructing the meson-exchange 
current is derived from the Lagrangian densities 

Inip 

C^u^A^) = -^-^e^p,,d''u;%x)d''u^'{x)r,'{x) , (15) 

where ect|3u^l is the Levi-Civita antisymmetric tensor with £0123 = +1- The vector meson- 
exchange current is then given by 

E {[r^^.(A^.)]i^/^a/3(A^.)rX.(A:.,/c;)^^M.(/c;)[r5^^.(^;)]2} , (le) 

V=p,U) 

where Dapik^) and Dfj_i,{k'^) stand for the (Feynman) propagators of the two exchanged 
vector mesons (either the p or u mesons as f = p or cu) with four-momentum and k'^, 
respectively. The vertices F involved in the above equation are self-explanatory. The A^A^u 
vertex Fj^^^(f = p, u) is taken consistently with those in the potential used for constructing 
the A^A^ T-matrix in Eq.(|). 

The coupling constant g^^ri' is determined from a systematic analysis of the radiative 
decay of pseudoscalar and vector mesons in conjunction with vector- meson dominance. This 



is done following Refs. |3ll , p!8| , with the aid of an effective Lagrangian with SU(3) flavor 
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symmetry and imposition of the OZI rule. The parameters of this model are the angle 
av{c(p), which measures the deviation from the vector (pseudoscalar) ideal mixing angle, 
and the coupling constant of the effective SU(3) Lagrangian. They are determined from a 
fit to radiative decay of pseudoscalar and vector mesons. The parameter values determined in 
this way in Ref. (model B), however, overpredict the measured radiative decay width of 
the 7]' meson |^2|. Therefore, we have readjusted slightly the value of the coupling constant of 



the SU(3) Lagrangian in order to reproduce better the measured width. We have ay = 3.8° 
and ap = —45°, as given by the quadratic mass formula, and the coupling constant of the 
effective SU(3) Lagrangian of G = 7 in units of 1/ \Jm„m'^^ where and stand for the 
mass of the two vector-mesons involved. The sign of the coupling constant G is consistent 
with the sign of the p7r7 and uji^^ coupling constants taken from an analysis of the pion 



photoproduction data in the ~ 1 GeV energy region ||33|. With these parameter values we 
obtain 

9ppri' = -G'sin(ap) = 4.95 

gujiur]' = —G (^V2 sin^ (ay) cos{ap) + cos'^{av) sin(ap)^ = 4.90 . (17) 
The VV7]' vertex {v = p,uj) in Eg. ([T6|) , where the exchanged vector mesons are both off 



their mass shells, is accompanied by a form factor. Following Ref. |18|, we assume the form 



^„,(e*-)^(|^)(^). (18, 

It is normalized to unity at = ml and k'^"^ = 0, consistent with the kinematics at which 
the value of the coupling constant gmr^i was determined. We adopt the cutoff parameter 



value of = 1.45 GeV as determined in Ref. |T8[ from the study of the u and (j) meson 
production in pp collisions. 

Another potential candidate for mesonic current is the cxT^T^'-exchange current, whose 
coupling constant may be estimated from the decay width of rj' into an rj and two pions. We 
take the Lagrangian densities 

(x) = -^^a{x)d,v{x)d^v'ix) 
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C...{x) = ^^a{x)d,^{x) ■ 9^7r(x) , (19) 

where cr(x) and ri{x) stand for the a and t] meson field, respectively, and m^, and m^' 
stand for the masses of the vr, rj and rj' meson. With the coupling constant of gan-n ~ 1-3, 



extracted from Ref. |24[ in conjunction with the NNa coupling constant used in the A^A^ 
potential V , and assuming that the measured decay width of r(?7' 77 + 27r°) = A2KeV 
is entirely due to the a meson intermediate state, we obtain an (upper) estimate of 
\9ar]ri'\ ~ 1.14. The sign of this coupling constant is not fixed. We shall consider both 
possibilities. 

The arjrj' current reads 

Jariri' = { [rATTv^] (g^) T^^^/ (g^, g^)i (gj) [Fiv (g^)] 2 } + (1^2) (20) 

in our previously-defined notation. The vertices T]\fNa- and Tn^^i are taken consistently with 
those in the NN potential V in Eq.(||). The crr/r/' vertex, ro.^^/(go., g,,), is provided with the 
same form factor given by Eg . ([ISD , except for the replacement m„ — > rrir,. 
The total mesonic current is then given by 

There are, of course, other possible mesonic currents, such as the uxprj'- and ^^r^'-exchange 
currents, that contribute to 1]' meson production in pp collisions. Their contributions can 
be estimated in a systematic way following Ref. |l^ and are found to be negligible. 



IV. NUMERICAL RESULTS AND DISCUSSION 

Once all the ingredients are specified, the total cross section for the reaction p + p ^ 
p + p + rj' can be calculated. We will first analyze each current separately, exploring the role 
of individual contributions and the uncertainties associated with them, and then combine 
them later in the section. 
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The A^* resonance current contribution to the total cross section as a function of excess 
energy Q is shown in Fig. ^. The excess energy is defined as Q = y^— where denotes 
the total center-of-mass energy of the system and = 2m at + m^/, its r^'-production 
threshold energy. Here we also display the recent experimental data from the COSY-11 
and SPESIII 0] collaborations. The mixing parameter A in Eq.(|TUp for the NN*tt vertices 
is fixed to be A = 0. Hereafter A^* stands for both the S'ii(1897) and Pii(1986) resonances 
unless otherwise indicated. The upper panel shows the results with A = in the NN*ri' 
vertices. We see that the cross section is largely dominated by the 5*11(1897) contribution 
(dash-dotted line). The Pii(1986) resonance leads to an energy dependence (dashed line) 
of the cross section steeper than that of the 5*11 (1897) resonance, although in the excess 
energy region considered its contribution is practically negligible. The solid line corresponds 
to the total contribution, which is close to the data. The lower panel corresponds to the 
choice of A = 1 at the NN*7]' vertices. Compared to the case of A = 0, we see that the 
mixing parameter has practically no influence on the 5*11(1897) contribution. This is due 
to the fact that this resonance is practically on its mass shell in the near-threshold energy 
region, so that both choices of the parameter A are equivalent. Although still very small 
compared to the 5*ii(1897) contribution, the Pii(1986) resonance contribution increases 
substantially in the low excess energy region compared to the A = case. This is because 
this resonance is off-shell and the ps coupling mixes efficiently the positive and negative 
energy states. The pv coupling (A = 0) suppresses this mixing. The interference with the 
5*11 (1897) contribution is now destructive. In this case the total contribution (solid curve) 
is slightly reduced compared to that in the upper panel. These results are consistent with 
the recent analysis of rj' photoproduction ||T6[. We recall that the coupling constants at the 



NN*ri' and NN*7[ vertices entering in the resonance current have been determined from the 
decay widths extracted in Ref. |T^ . 



Although we find the consistency between our result above for pp —>■ pprj' and the rj' 
photoproduction analysis of Ref. |]IB[, care must be taken in drawing any premature con- 



clusions about the basic rj' production mechanisms in these reactions. As mentioned in 
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subsection III.B, given the assumptions made in Ref. in the extraction of the partial 



decay widths, the corresponding couphng constants, gNN*ri' and gNN*n, are subject to con- 
siderable uncertainties. In particular, the resonance dominance assumption does not allow 
any other possible production mechanism. As we shall show below, the vector-meson ex- 
change current yields also a cross section comparable to the data. This indicates that the 
resonance current contribution shown in Fig.|l] is probably overestimated. Another source 
of uncertainty in the resonance current arises from the off-shell form factors introduced at 
the NN*iT vertex. Although the dominant 5*11(1897) resonance contribution is insensitive 
to the form factor at the NNr]' vertex because the resonance is nearly on its mass shell 
near the threshold energy, the corresponding exchanged pion (see Fig.|l]) is far off-shell. This 
makes the resonance contribution very sensitive to the form factor involving this pion. In 
the present calculation we have employed the same monopole form factor used at the NNtt 
vertex entering in the NN potential V in Eq.(0) with the cutoff parameter of A^ = 0.8 GeV. 
One might then argue that, within the resonance dominance assumption, the jp —>■ pr]' and 
pp pprj' processes can be described consistently — provided that one uses a soft form fac- 
tor at the NN*iT vertex. Whatever the arguments, it is clear that further investigation is 
required in order to determine better the relevant coupling constants before a more definite 
conclusion can be drawn about the role of the nucleon resonance current in the pp pprj' 
reaction. Of course one should also keep in mind the possibility that the strong peak in the 
cross section observed in rj' photo-production close to threshold may be due to a cusp effect 
p9| and not due to a resonance. 

We turn now to investigate the mesonic current contribution, which is displayed in Fig. |^. 
Among the exchange currents considered, the pprj' exchange gives rise to the dominant 
contribution (dashed line). The ujurj'- and ar^r^'-exchange current contributions (dotted and 
dash-dotted line, respectively) are of similar magnitude. However they are smaller by about 
a factor of five compared to the ppr/'-exchange contribution. The slightly different energy 
dependence of the cross section resulting from the ppr]'- and cucJT^'-exchange currents is due 
to the tensor coupling in the NN p vertex. This coupling is absent in the NNuo vertex 
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used here. As mentioned in subsection III.C, the sign of the (T?7?7'-exchange current relative 
to f f r^'-exchange currents [v = p, u) is not fixed; therefore we consider both possibihties. 
Assuming the negative sign for garjr]', the interference with the f f r/'-exchange current is 
destructive, yielding the total contribution represented by the solid line. If g^r^r]' is positive, 
then the interference is constructive and gives the total contribution shown by the dotted 
line. Although this latter choice leads to an overprediction of the data, we cannot exclude it 
at this stage because, as we shall show below, the nucleonic current can interfere destructively 
with the mesonic current so that the combined contribution may be smaller than that of 
mesonic current alone. Moreover, the off-shell form factor at the r^'-production vertices in 
the mesonic current is not well-determined at present. As mentioned in subsection III.C, we 
employ the same form factors at the npu and 7rp0 vertices as in the study of u- and 0-meson 



production in pp collisions [|18|, which are themselves subject to considerable uncertainties. 
Indeed, a reduction of about 10% in the value of the cutoff parameter can bring the short- 
dashed curve in Fig. ^ onto the data. Therefore neither sign of (7^777)' can be excluded from 
the present analysis. In this regard, measurements of cross sections at higher excess energies 
may be useful in deciding between the two choices of the sign in the arji]' coupling. 

In the upper panel of Fig. ^, the nucleonic current contribution (dashed curve) is shown 
with pv coupling at the NNrj' vertex. The contribution from the mesonic current with the 
choice of positive sign for g^j^^i' (dashed curve in Fig. ^ is also shown here as the dash-dotted 
curve. As we can see, the nucleonic current is about two orders of magnitude smaller than 
the mesonic current. Their interference is constructive and results in the total contribution 
shown by the solid curve, which overpredicts the data. The lower panel of Fig. ^ shows the 
same results, except that ps coupling is used at the NNrj' vertex in the nucleonic current 
contribution. Compared to the case of pv coupling, we observe a considerable enhancement of 
the nucleonic current. This is due to the strong admixture of the positive and negative energy 
nucleon states provided by the ps coupling; with pv coupling this mixing is suppressed. 
Furthermore the interference with the mesonic current is now destructive, resulting in a total 
contribution (solid line) that is substantially smaller than the mesonic current contribution 
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alone. Fig. ^ shows the same results as in Fig. |, except that here the mesonic current 
contribution corresponds to the choice of negative sign for g^r^f^i (solid curve in Fig. We 
see a similar feature to that observed in Fig. yielding net cross sections that underpredict 
the data. In particular, we see a strong destructive interference between the mesonic and 
nucleonic currents in the case of the ps coupling at the NNrj' vertex (lower panel). We 
observe that the effect of the nucleonic current may be overestimated here because we have 
used a rather large NNrj' coupling constant, as discussed in subsection III. A. In any case, 
given the fact that the nucleonic current is small compared to other currents, a quantitative 
determination of the NNrj' coupling constant using the pp — > pprj' reaction is unlikely to be 
possible — at least not before the other dominant currents are better understood. 

The above considerations expose the role of individual currents as well as the current 
uncertainties associated with them in the description of the pp pprj' reaction. We combine 
now all the currents. In the following we consider all the possibilities for the mesonic-plus- 
nucleonic current contribution as shown in Figs. ^ and |^. For each scenario considered we 
adjust the coupling constants gNN*ri' and gNN'-K in the resonance current so as to reproduce 
the data, taking A = 1 in the NN*t]' vertices. The results are shown in Fig. ^ For the 
case of the mesonic-plus-nucleonic current shown in the upper panel of Fig. ^ the corre- 
sponding resonance contribution is set to zero because the mesonic-plus-nucleonic current 
already overpredicts the data and the resonance current can only enhance the cross sec- 
tion further. We recall that the interference between the resonance and other currents is 
practically negligible. Therefore we have here adjusted the ps-pv mixing parameter in the 
NNrj' vertex to A = 0.7 in order to reproduce the data, and the corresponding result is 
shown in the solid curve of Fig. ^ For the choice of the mesonic-plus-nucleonic current as 
shown in the lower panel of Fig. ^ the resonance current is required in order to reproduce the 
data. We have adjusted this current contribution by multiplying the product of the coupling 
constants gNN*ri'gNN*TT for both A^* = 5*11(1897) and Pii(1986) resonances by an arbitrary 
reduction factor of 1/2. The result is shown by the dash-dotted line. For the scenario of 
the mesonic-plus-nucleonic current shown in the upper panel of Fig. |, the reduction factor 
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of the product gNN'q'QNN'-K required in the resonance current is 0.625. The corresponding 
result is shown by the dotted curve. Finally, for the scenario of the lower panel in Fig. 
no reduction of the resonance current is required in order to reproduce the data, since this 
corresponds to an extreme case where the destructive interference between the nucleonic 
and mesonic currents leads to a mesonic-plus-nucleonic current contribution that is nearly 
negligible. The short-dashed curve is the corresponding result. Fig. |^ reveals the different 
energy dependences arising from each scenario considered. We see that if one wants to learn 
about the production mechanism from the total cross section then one has to go to much 
higher energy region than that covered by the presently available data. 

At this point we note that, in view of the latter scenario discussed in Fig.|| above, it 
is not surprising that the approach of Ref. §| describes the pp — > pprj' data. In Ref. 
the TiN Tj'N amplitude, extracted from the measured cross section, has been used as 
the production mechanism. The vrA^ rj'N amplitude obtained in this way must include 
the nucleon resonance contribution that we have considered in this work, which leads to a 
strong vrA^ — > r]'N transtion. We emphasize, however, that one has to be careful in drawing 
any premature conclusion here. First, as discussed before, our nucleon resonance current is 
subject to considerable uncertainties. Second, although the calculation of Ref. |^ is based on 
the measured ttN — » Nrj' cross section, it requires an introduction of an off-shell form factor 
at the NNn vertex in order to account for the off-shellness of the exchanged pion. As we 
have pointed out, the exchanged pion is far off-shell in the present reaction and consequently 
the resulting cross sections are very sensitive to the corresponding form factor at the NNtt 
vertex. We note that in Ref. a monopole form factor with the cutoff parameter of 
Att = 1.3 GeV, as used in the Bonn A^A^ potential, has been employed. There are, however, 
a number of indications that the NNtc form factor is much softer (Ayr ~ 0.8 GeV) than that 
used in the Bonn potential p5|j36|] . This would reduced the contribution of the one-pion 



exchange to the pp — > pprj' process substantially. The above consideration gives an idea of 
the kind of uncertanties involved in the theoretical predictions. 

We now focus our attention on the excess energy region close to threshold where the data 
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exist. Fig. ^ shows the results of Fig. ^ in the range of excess energy up to Q = 10 MeV. 
Here, only one total current contribution (solid curve in Fig. |^) is shown since all four sce- 
narios considered lead to almost indistinguishable results at the scale used in the limited 
energy domain considered in the figure. We note that the energy dependence of the cross 
section is not reproduced. In particular, after adjusting the ps-pv mixing parameter or the 
reduction factor of gNN'ri'dNN'-K in order to reproduce the higher energy data points, as has 
been done in Fig. ||, we are not able to describe the lowest two data points. We observe that 
the data have considerable uncertainties in the excess energy, presumably due to background 
subtraction that becomes more difficult for energies very close to the threshold. Nevertheless 
the results indicate the need for some mechanism(s) that introduces an additional energy 
dependence of the total cross section very close to threshold. The energy dependence of the 
total cross section near threshold is basically determined by the final state A^A^ interaction 
T^ . This is illustrated by the dotted line in Fig. |^, which corresponds to the result calcu- 



lated taking into account only the three-body phase space-plus-final state interaction; it has 
been normalized to the solid curve in the low excess energy region. We see that the pro- 
duction currents start to modify the energy dependence given by the phase space-plus-final 
state interaction only for energies above (5 = 4 — 5 MeV. Therefore the predicted energy 
dependence in the near- threshold region is unlikely to be modified by the basic production 
mechanism(s). In fact, all the individual currents considered in this work lead practically to 
the same energy dependence in the restricted energy range considered in Fig. 0. 

The fact that the energy dependence given by the three-body phase space-plus-final state 
interaction does not reproduce the observed energy dependence in pp — > ppr]' process has 
been a subject of attention since the publication of the COSY-11 data P,p|,p!T| . In Fig. | we 



show the measured cross section data for tt°, rj and rj' production in pp collisions together 
with the corresponding energy dependence predicted by the phase space-plus-final state 
interaction only (solid curves). Apart from the fact that in rj' production the deviation 
starts a little lower in excess energy than in tt° production, we see no indication for any 
peculiar feature in rj' production. 



The most trivial source of influence on the energy dependence of the cross section close 



to threshold is the Coulomb force [0, which has not been considered so far in our calcula- 



tions. We follow the Gell-Mann and Goldberger two-potential formalism to include the 
Coulomb force in the flnal state interaction. The perturbed wave function is calculated ex- 
actly by solving the Schrodinger equation in coordinate space with the nuclear-plus-Coulomb 
potential. For the nuclear potential we use the Paris potential . The half-off-shell (pure) 



Coulomb T-matrix is obtained in analytic form [^. We note that care must be taken in 



order to use consistently a NN T-matrix constructed within a non-relativistic formulation- 



such as the Paris T-matrix — in a relativistic approach. We have followed Ref. in order 
to use the Paris T-matrix in the present relativistic approach. Also, in the following, the 
production current has been kept the same as that used in conjunction with the Bonn B 
potential, except for a readjustment of the reduction factor of the product gNN*r]'gNN*n in 
the nucleon resonance current that was needed to reproduce the data. We are aware that by 
using the Paris potential instead of the Bonn B potential we may loose the consistency be- 
tween the production current and the flnal state interaction, however the energy dependence 
of the cross section is practically the same — when the Coulomb force is switched off — as that 
given by the Bonn B potential over a wide range of excess energy. The result of including 
the Coulomb effect as described above is shown in Fig. || by the solid line. The dotted 
line corresponds to the result when the Coulomb force is switched off. As can be seen, the 
Coulomb effect is considerable; in particular, it reduces the cross section by almost 40% at 
Q = 1.5 MeV. 

The picture that emerges from the above considerations is that there is no obvious indica- 
tion of a need for other mechanisms than those already considered here in order adequately 
to describe the existing data. A close comparison of the COSY-11 and SPESIII |^ data, 
however, seems to indicate different trends in the energy dependence. In this connection, 
new measurements at higher energies, and even remeasurement of some of the data near 
Q ~ 4 MeV, are extremely important for resolution of this issue. The new data at higher 
energies will certainly impose more stringent constraints on the parameters of our model. 
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They would also tell us whether a proper description of the pp pprj' requires the inclusion 
of the prj' final state interaction, as has been argued recently [|11|] . The prj' final state inter- 



action as discussed in Ref. |]TT| can be viewed as the higher order terms in the production 
current treated in the present work. In fact, the ?7'-emission current defined in the previous 
section (cf. Eqs. (p|,P|,pTD ) is, in part, just the Born term of a more general current that can 



be obtained by using the T-matrix amplitude for the M + p ^ rj' + p transition, where M 
denotes any meson of interest. This can easily be seen if we disregard the nucleon labeled 2 
in Fig. the current then becomes nothing other than the Born term of the M + p rj' +p 
transition amplitude. The T-matrix amplitude for the M + p — >■ rj' + p process may be 



separated into the so-called pole and non-pole terms according to Pearce and Afnan 
By using the physical baryon mass and physical NBr]' {B = N, N*) and NNM vertices 
in the baryonic current, the pole term of the M + p ^ rj' + p T-matrix amplitude is fully 
accounted for in the present work. What is taken in the Born approximation is, therefore, 
the non-pole part of the T-matrix only. 

V. SUMMARY 



We have investigated the reaction pp pprj' within a relativistic meson-exchange model 
of hadronic interactions. We find that the 5*11(1897) resonance as determined in a recent 
multipole analysis of rj' photoproduction |jl6| gives rise to a contribution that reproduces the 



measured cross section and allows for no other production mechanisms. This is consistent 



with the analysis of Ref. Il^] , where the nucleon resonance parameters have been extracted 
under the resonance dominance assumption. Whether this is already the entire physics of 
the production mechanism is an open question that requires further investigation. Indeed, in 
Ref. |T^, in addition to the assumption of resonance dominance, only the prj' and pvr" decay 
channels have been explicitly considered; the pvr" decay channel accounts for all possible 
decay channels other than the prj' . The coupling constants NN*ri' and NN*tt obtained from 
the decay widths extracted under these assumptions are, therefore, subject to considerable 
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uncertainties which, in turn, lead to corresponding uncertainties in the resonance current 
contribution to the cross section. A more accurate determination of the relevant resonance 
parameters is called for before a more definite conclusion can be drawn about the role of 
the N* resonances in the pp — > pprj' reaction. We should also keep in mind that the strong 
peak in the cross section close to threshold in rj' photoproduction may be due to a cusp 
effect [|^, and not to the resonance assumed in Ref. [|l6l. All these issues associated with 



the resonance require further careful investigation. 

The mesonic current gives rise to cross sections which are comparable to the measured 
values. It is dominated by the ppr/'-exchange current. Depending on the choice of the 
sign of the arjr]' coupling constant, the calculated cross sections overpredict or underpredict 
the data. However neither choice can be excluded at present — mainly due to uncertainties 
associated with the baryonic current. Data in higher energy regions than currently available 
may provide the necessary constraint to fix the sign of this coupling, as the different signs 
lead to different energy dependences. 

Since the nucleonic current is found to be small, it is unlikely that the pp — > ppt]' reaction 
can be used quantitatively to determine the NNrj' coupling constant — at least not before 
the other dominant currents are under better control. 

As expected from earlier calculations [|l3l, the Coulomb force in the final state inter- 



action is found to play a crucial role in explaining the observed energy dependence of the 
cross section close to threshold. The basic production mechanisms start to influence the 
energy dependence of the cross section only for excess energies above ~ 5 MeV . For lower 
energies the energy dependence is determined by the three-body phase space-plus-final state 



interaction [0. New data at higher energies will certainly impose constraints on the pro- 
duction mechanisms. In addition to more exclusive observables than the total cross section, 
a combined analysis of production processes using both electromagnetic and hadronic probes 
should help to disentangle the different production mechanisms. 
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FIG. 1. r/' meson production currents included in the present study. Upper row: nu- 
cleonic current Jnwo ^ = 7!", ?7, p, w, cr, Oo- Middle row: nucleon resonance current Jres, 
N* = S'ii(1897), Pii(1986). Lower row: mesonic current Jmes, v = p,LO. 
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FIG. 2. A^* resonance current contribution to the total cross section for the reaction pp ppij' 
as a function of excess energy Q. The upper panel corresponds to the results with the choice A = in 
Eq.(p^) at the NN*r]' vertices, the lower panel to those with the choice A = 1. The dashed-dotted 
curves correspond to the 5ii(1897) resonance contribution, the dashed curves to the Pii(1986) 
resonance contribution. The solid curves are the total contributions. The experimental data are 
from Refs. IQ. 
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FIG. 3. Mesonic current contribution to the total cross section for the reaction pp ppij' 
as a function of excess energy Q. The dash-dotted curve correspond to the o"???7'-exchange and 
the dotted curve to the lolot]' -exchange current contribution. The dashed curve represents the 
ppr/'-exchange current contribution. The sohd curve is the total contribution corresponding to 
choosing the negative sign for the ar]?]' coupling constant, the short-dashed curve to the positive 
sign. The experimental data are from Refs. §,0]. 
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FIG. 4. Nucleonic-plus-mesonic current contributions to the total cross section for the reaction 
pp — > ppr]' as a function of excess energy Q. The upper panel corresponds to the results with 
the pv coupling (A = 0) at the NNt]' vertex in the nucleonic current, the lower panel to those 
with ps coupling (A = 1). The dashed-dotted curve in both panels represents the total mesonic 
current contribution shown by the dashed curve in Fig. |3|. The nucleonic current contributions are 
represented by the dashed curves and the corresponding total contributions by the solid curves. 
The experimental data are from Refs. [^,@]. 
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FIG. 5. Same as Fig. ^, except for the total mesonic current contribution, which corresponds 
to the result represented by the solid curve in Fig. ^ 
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FIG. 6. The total cross section for the reaction pp — > ppr]' as a function of excess energy 
Q. The sohd hne represents the mesonic-plus-nucleonic current contribution corresponding to 
the upper panel of Fig. ^, except for the value of the ps-pv mixing parameter at the NN-q' ver- 
tex, which is taken to be A = 0.7. The dotted line represents the result corresponding to the 
mesonic-plus-nucleonic current shown in the lower panel of Fig. ^; it includes also the nucleon reso- 
nance current, whose product of the coupling constants gNN*r)'9NN*-K is multiplied by an arbitrary 
reduction factor of 1/2. The dash-dotted and dashed curves are the same as the dotted curve, 
except that they correspond to the mesonic-plus-nucleonic current shown in the upper and lower 
panels of Fig. |5|, respectively. The corresponding reduction factors of the product gNN*Ti'9NN*w are 
0.625 and 1. The experimental data are from Refs. [P,|7|. 
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FIG. 7. Same as Fig. ^ for excess energy below 10 MeV. Only one total contribution (solid 
curve) is shown since all the curves shown in Fig. ^ are almost indistinguishable at the scale 
used here. The dotted line corresponds to the result with only the phase space-plus-final state 
interaction. 
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FIG. 8. The total cross section for tt"-, rj- and r/'-meson production in pp collisions as a function 
of excess energy Q. The curves represent the corresponding energy dependence given by the phase 
space-plus-final state interaction with no Coulomb force. In the upper panel the experimental data 
are from Refs. |42-Q; in the middle panel from Refs. ||^,^,4(:]; in the lower panel from Refs. ||6|,^. 
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FIG. 9. The total cross section for the reaction pp pp-q' as a function of excess energy Q. 
The sohd curve corresponds to the total current contribution including the Coulomb force in the 
final state interaction. The dotted line corresponds to that when the Coulomb force is switched 
off. Here the Paris T-matrix has been used as the final state interaction. The experimental data 
are from Refs. ||6| (circle) and (triangle). 
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